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The Majorana resonant level model (MRLM) describes the universality class of the two- 
channel/terminal Kondo model at the Toulouse point as well as of a resonant level between two 
half-infinite Tomonaga-Luttinger liquids. We analyze the time evolution of the electric current and 
of the population function after an instantaneous switching on of the tunneling coupling. We find 
that the only timescale, which governs the relaxation of the initial dot preparation is the inverse 
contact transparency F, whatever the dot offset energy A, applied bias voltage or temperature. The 
voltage alone determines the superimposed oscillatory behavior of the observables for weak detun- 
ing |A| < r/2. In the opposite case of strong detuning |A| > F/2 a beating pattern emerges. For 
the current the finite temperature plays the similar role as the hybridization. The dot population 
function dynamics approaches that of a resonant (A = 0) setup upon increasing the voltage or/and 
temperature. 

PACS numbers: 72.10.Fk, 73.63.-b, 71.10.Pm 
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During the last decades we witnessed a soaring inter- 
est in the nonequilibrium properties of quantum impu- 
rity models. It has been motivated by the advances in 
manufacturing of semiconductor based nanostructures, 
in which only a small number of electronic levels par- 
ticipate in the charge and spin transport. The exper- 
imentalists not only succeeded in observing the Kondo 
effect in such structures, but also were able to demon- 
strate phonon-assisted charge transport through single 
molecules, see e. g. 0. From the theoretical point of 
view such zero-dimensional systems can be seen as real- 
izations of quantum impurity models with or without in- 
ternal degrees of freedom. Typical examples are the An- 
derson model, Tomonaga-Luttinger liquids (TLLs) and 
fractional quantum Hall edge states with impurities, res- 
onant levels coupled to a Holstein phonon etc. Many of 
their equilibrium properties can be calculated with the 
help of by now quite substantial number of techniques. 
On the contrary, in the nonequilibrium one is especially 
missing the extremely powerful Bethe ansatz and integra- 
bility methods. Although their adaptation to setups far 
from equilibrium has already begun, a number of open 
issues are yet to be settled. [2| 

Fortunately, at least in the universality class of Kondo 
models there are some exactly solvable points which oc- 
cur at nontrivial constellations of system parameters: the 
Toulouse limits. Those turn out to be the only fully fea- 
tured testing grounds (apart from numerical approaches, 
of course), which we have at our disposal in the mo- 
ment. In its canonical incarnation it is a special point 
in the parameter space where a mapping of the con- 
ventional Kondo model onto a resonant level model is 
possible. 0, Q] For transport phenomena a two terminal 
topology is necessary. [B[ The original Hamiltonian is then 
mapped onto that of a Majorana resonant level model 
(MRLM). It turns out to be isomorph to the interact- 
ing resonant level in the TLL.Q Hence, a solution of 



MRLM provides us with predictions to a considerably 
broader class of situations, than the original Kondo quan- 
tum dot. Thus far this kind of exact solvability has been 
taken advantage of in a number of contributions. 0, @, @] 
However, only a fraction of them went beyond the steady 
state calculations, ljl U| The quantities of special inter- 
est are the transport current evolution, several aspects of 
which has bee discussed in Ref. [1(3], and the magnetiza- 
tion of the impurity in the Kondo case, or the population 
of the resonant level in the TLL resonant level setup. To 
the best of our knowledge, the latter observables have 
not yet been addressed. We would like to close this gap 
and present the full solution in the case of instantaneous 
switching on of the tunneling hybridization. 

We first summarize the MRLM Hamiltonian in its most 
general form. It features local Majorana fermions a, b as 
well as four different Majorana fermionic fields. They 
are organized into pairs describing a principal channel 
rj(x),^(x) and a, flavor channel rjf(x),£f(x): 

H = H [t,r],Zf,fif]-i[Aab + J-btf{0) + J + ofif(0) 
+ 7+ H(0)+7-ar 7 (0)] . (1) 

Hq is responsible for the dynamics of the free (decoupled 
from the local Majoranas) fields, 



H 



dx 



r]f(x)d x r)f(x) + Zf(x)d x £ f (x) (2) 



ri(x)d x rj(x) + £(x)d x Z(x) + V^(x)r](x) 



A, J± and j± arc different constant couplings whereas V 
is in a strict sense not a coupling but a chemical potential 
and is related to the bias voltage applied in the systems 
parental to MRLM. This is the reason why we employ 
nonequilibrium diagrammatics in order to calculate the 
observables of interest. These are the current 



J=(i/2)[ 7 _ (o£(0)> + 7+ <6 »?(())>] 



(3) 
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and the dot occupation probability, 
n d = (l + i(a6))/2. 



(4) 



The origin of the MRLM is twofold: it emerges in the 
(i) two terminal Kondo model at the Toulouse point, 
and (ii) interacting resonant level between two TLLs. 
In the situation (i) one starts with the conventional 
Kondo Hamiltonian (we set h = Vp = e = fcg = 1), 
H = H + Hj + Hm + Hy, where, with ip a <T being the 
electron field operators in the R,L (right/left) terminals, 

i?0 = ! ^ E / dx ^ta( x )dx^aa{x) , 
a=R,L cr=T,| ^ 

#7 = E E 

a,0=R,L v=x,y,z 

Hy = (y/2)E / dx(^ La -^ R J JRa ), 
H M = -»B9iHT z = -At z . (5) 

Here t v=x ' v ' z are the Pauli matrices for the impurity spin 
and (a, (3 — R, L; a =|, j; er^, are the components of the 
vth Pauli matrix) 

are the generalized electron spin densities in (or across) 
the leads biased by a finite V. The last term in Eq. ([5|) 
stands for the magnetic field, A = [IbQiH. Following 
Ref. 0, we assume = jf = jf , J^ L = J RR = J z 
and = Jf L = 0. After bosonization, Emcry- 

Kivelson rotation, and refermionization (see details in [7] 
or [ijj]) and setting J z = 2it, one obtains the Toulouse 
point Hamiltonian (p}, where J± = ( Jf L ± J RR ) /y/2irao, 
7+ = J± = J RL / y/2na n , 7_ = (ao is the lattice con- 
stant of the underlying lattice model) and a and b being 
local Majorana operators originating from the impurity 
spin r = (t x ,t v ,t z ): t x — a, t v = b. The fields rjf and 
£f in the spin-flavor sector are equilibrium (real) Ma- 
jorana fields, whereas 77 and £ in the charge-flavor sec- 
tor are biased by the transport voltage V. The current 
through the system is then given by (O and the impurity 
magnetization by 



( T z ) = -i(r x t v ) = i (a b) = 2n d - 1 



(6) 



The setup (ii) is the spinless TLL resonant level model 
with the Hamiltonian H = Hk + H t + Hj, where Hk is 
the kinetic part, Hk — A d) d+^ i=R L Ho[ipi], describing 
the electrons in the leads Ho[ipi], and the resonant level 
with energy A, the corresponding electron operators be- 
ing d\d. The dot can be populated from either of the two 
leads (i = R, L) via electron tunneling with amplitudes 
7i, H t = 2i7i[rf^i(0) + h.c.]. Hj describes the electro- 
static Coulomb interaction with the strength U between 



the leads and the dot, Hi = Ud^d £V $(0)^(0). The 
contacting electrodes are one-dimensional half-infinite 
electron systems. We model them by chiral fermions liv- 
ing in an infinite system. Q In the bosonic representation 
Ho{ipi] are diagonal even in presence of interactions (for 
a recent review see e.g. 0] ; we set the renormalized Fermi 
velocity v = vp/g = 1, the bare velocity being vp): 



HobPi] = ^ J dx Px4(a 



(7) 



Here the phase fields <j>i(x) describe the slow varying 
spatial component of the electron density (plasmons), 
tpl(x)tpi(x) — d x 4>i{x)/2-ny/g. The electron field operator 
at the boundary is given by ipi(Q) = e i( ^ )/ ^ g ~/\/2ira . 
Here g is the conventional TLL parameter related to 
the bare interaction strength Utll via g = (1 + 
^TLL/TTViO^^-lilH] I n the chiral formulation the bias 
voltage amounts to a difference in the densities of the 
incoming particles in both channels far away from the 
constriction. 14| The current is then proportional to the 



difference between the densities of incoming and outgoing 
particles within each channel. 

At generic g 7^ 1 the problem cannot be solved exactly 
though. However, at g — 1/2 after a transformation of 
d* and d operators to the spin representation of the form 
t x = (dt + d), t v = -i{d^ - d), and r z = 2dU - 1 
one immediately observes that the [/-term is analogous 
to the T z -spin density coupling in the Kondo problem. 
Indeed, a mapping onto the MRLM can be performed in 
the same way as in the Kondo case and one again ob- 
tains Hamiltonian ([I]) with y± = 71, ± 7_r and J± =0.0 
As before we are interested in the transport current ([3]) 
and dot population (|4|). The time derivative of rid is the 
displacement current in the system -^jjgp(i) = drid(t)/dt. 
The currents through the individual contacts Il.r can 
then be conveniently evaluated using the relations [lij 
lL,R.{t) — I(t) ± -^dispW/2 - The effect of the asymmetry 
7_ 7^ is similar to that of a finite detuning A ^ 0.fl2| 
that is why in order to make the calculations and re- 
sults more lucid we concentrate on the latter case and 
set 7_ = 0. Then the current operators are identical in 
both (i) and (ii). 

In general the calculation of the observables is sim- 
plified by their reduction to a related nonequilibrium 
Green's function (GF). In the case of the transport cur- 
rent it is G bri (t,t') = -i(T c b(t)r)(tf)), where T c is the 
time ordering along the Keldysh contour C.[l6| Then the 
time dependent current is given by I(t) = —(7/2) Gb v (t+ 
+ ,t), where by abuse of notation Gb v {t,t') denotes the 
time ordered GF and 7 = 7+. From now on we assume a 
steplikc switching on of the tunneling when j(t) = 7 0(f), 
where 0(t) is the Heaviside function. By an expansion 
in 7 and resummation of the series one can show, that 



■A 



the following reduction is valid: [if 



I(t)=ij I dt'[D bb (t,t')g (ri (t' -t) 
-D<(t,t')gf v (t'-t) 



(8) 



Dbb(t,t') = —i(Tcb(t)b(t')) is the exact homogeneous b- 
Majorana GF. g^(t,t r ) = —i(Tc£,(t)r](t')) are the zero 
order in tunneling GFs. Remarkably, all its Keldysh com- 
ponents are equal and given by g^ v (^>) — — n ii)/2, 
where tir.j, are the Fermi distributions in the right (left) 
electrode. [12| In the Kondo setup we model the electrodes 
by wide flat band with constant density of states po. In 
the resonant level TLL setup po is equal to the energy- 
independent prefactor in the density of states in vicinity 
of the Fermi edge. Here the applied voltage must be 
doubled. [T3| Then for the current through the system we 
obtain: 



rr 



sinh[7r2 (t — t)\ 



where V = po7 2 /2- Thus everything is determined by the 
retarded Dj/ b only. In the time domain we obtain: 

D b %{t,t>) = Df b )R {t-t>) (10) 

/•OO 

Jo 

where the self-energy is due to the tunneling only and is 
up to a prefactor identical to the homogeneous 77-GF. It 
is most compact in the energy representation: 



-ir, S < (lu) — — ir(n^ + nn) , 
-iT(n L + n R - 2) . (11) 



To make progress we use 



D 



(0)R 
bb 



(t, t') = -iO(t - t') cos[A(i - t')] , (12) 



for the GF in the absence of coupling, which is valid at 
J± = 0.[l9( It can be shown, that the solution of (flQ|) 
is translationally invariant in time and has the structure 
t') = Kit - t>) = -»0(t - t') f(t - t'), where 



/(t) = /(0)(t)-r / drf(°\t-r)f(r), (13) 



with /(°)(t) = cos(Ai). This is a Volterra integral equa- 
tion of the second kind solvable by the Laplace transfor- 
mation. We observe however, that the equation for the 
retarded dot GF in the static case without switching ef- 
fects has exactly the same form.[20( The explanation for 
that is the functional similarity of our switching method 
and of the natural steplike time dependence of the re- 
tarded GF. The solution in the energy domain is known 
to be given by 



j 2 - A 2 + iuT 



(14) 




— v=ior, a=o 
■ ■ ■ ■ v=2r, A=0 

— - v=ior, A=4r 

— v=ior, A=ior 
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FIG. 1: Full current through the constriction at zero temper- 
ature for different voltages V and dot offset energies A. The 
current is measured in units of Go T/2, where Go = 2e 2 /h is 
the conductance quantum. 



Transformed back it yields 

-rt/2 

/(<) = 2 ^ [2n cosh (ifi) - T sinh (ifi)] , (15) 



with Q. = ,/(r/2) 2 - A 2 for weak detuning |A| < T/2. 
In case of strong detuning, when |A| > T/2 the function 
is found by an analytic continuation. Knowing that the 
retarded GF is indeed translationally invariant in time 
allows for further simplification of ((9|): 



m = ™ I* sin(W) 
W 2 ] Jy ' sinh(TrTr) 



(16) 



In general the current shows up two constituents: the 
transient one and the time- independent static one, I(t) = 
^stat + ^trans W- For the weak detuning we obtain[6j,[l0| 



4stat 



— Im5> +P r/2) 

p=± 

1 iV + pn + T/2" 

2 + 2VT 



(17) 



where ^ denotes the digamma function. The transient 
part is given by 



'trans 



(t) 



TT T ^ (n -pT/2) e dv+ P n-r/2-irT)t 

— — Im > 

2VL ^1 

p=± 



x 2 Fj ( 1 

3 
2 



1 



iV +pfi-r/2-7rT 
iV+pQ-T/2 



iV - 



2-kT 

■pa -T/2 -2nTt 

1 " 



2irT 



(18) 



where 2F1 denotes the hypergeometric function. For 
weak detuning the transient current oscillates with the 
frequency cx T^.[l^ This situation changes in case of 
strong detuning when a becomes imaginary and enters 
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FIG. 2: Population of the dot at zero temperature for dif- 
ferent voltages V and dot offset energies A. The dotted line 
corresponds to the resonant case A = 0. 



the above equations in the same way as the voltage. Then 
there are two different frequencies oc |y±f2| and a beating 
pattern emerges, s. Fig. [TJ Of course these features can 
be observed only when the oscillation period is smaller 
than the competing time scale oc T + ttT, which governs 
the overall current relaxation. 

It turns out, that for the calculation of the dot popula- 
tion it is more convenient to work with the lesser Keldysh 
GF D< b (t, t') = -i(a(t) b(t')). The calculation of this GF 
is accomplished using the Dyson equation in the time do- 
main (multiplication corresponds to time integrations): 



Dab 



D 



(0) 



where the self-energy is pro- 



portional to the unperturbed GF for the £-Majoranas 
£ = Tg^. After the Keldysh disentanglement we obtain 



d: 



D 



(0)< 



D 



^ u bb 



(19) 



, n (0)R y< D A , n (0)< y A n A 
+1J ab ^ U bb + U ab ^ U bb 



Thus the calculation of the inhomogeneous D a &-GF is 
now reduced to the calculation of the homogeneous D bb . 
The necessary zero order GFs are 

DTit) = "(«/2) e lKAt , D^ R (t) = -*0(t) Bin(At) , 

where k = ±1 encodes the initially populated/empty 
dot level. For the retarded component of the homo- 



geneous D bb one obtains the equation = D 



(0)R 
'bb 



u bb ^ u bb- 



(0)R . 
bb 

The resonant case A = is especially simple. Here 
the bare retarded inhomogeneous GF vanishes and one 
obtains 



nd(t) 



1 



ne 



-rt\ 



/2- 



(20) 



This is a very remarkable result since the relaxation of the 
dot occupation/magnetization is independent on both 
temperature and applied voltage and is only governed 
by the timescale 1/r. 



In the off-resonant case the calculations are more com- 
plex since we need the full homogeneous D^. It is given 

by0 



D<=(l + D' b l^)D^<(l 
Putting the result into (|T9|) we obtain for weak detuning 



^ A D£ b ) 



rt R v< n A 

u b b u u bb ■ 



, . 1 , rt , TA f duj 

nAt) = - (l + ne~ vt ) / — — — 

dy ' 2 V ' 20 2 J 2tt (A 2 



2p2 



i 



,-rt 



)-e 



-rt/2 



[{uo 2 + A 2 ) sinh(fit) sin(wt) 
+2u>n cosh(fM) cos(wt)]} . 



The most striking feature of this result is that the re- 
laxation behavior of the information about the initial 
preparation is exactly the same as in the resonant case 
and independent on either the temperature or the applied 
voltage. Thus the rate at which the system 'forgets' its 
initial preparation does not depend on these parameters. 
At zero temperature the effect of finite voltage is to fix the 
upper energy integration boundary. Since the integrand 
is an odd function of energy the whole correction due to 
finite A vanishes towards larger V, so that rid approaches 
that of a system at resonance, see Fig. [2j The effect of 
finite temperature is very similar: due to smearing off of 
the Fermi edge the contribution of the energy integra- 
tion decreases so that the relative weight of the resonant 
contribution (|20p increases and becomes more and more 
dominant. For finite detuning rid oscillates with the fre- 
quency oc A. As soon as the applied voltage becomes 
nonzero a beating pattern emerges just as for the time 
dependence of the transport current. 
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